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Abstract 

A new Diagonally Inverted LU Implicit scheme is de- 
veloped within the framework of the multigrid method for 
the three-dimensional unsteady Euler equations. The ma- 
trix systems that are to be inverted in the LU scheme are 
treated by local diagonalizing transformations that decou- 
ple them into systems of scalar equations. Unlike the Diag- 
onalized ADI method, the time accuracy of the LU scheme 
is not reduced since the diagonalizing procedure does not 
destroy time conservation. Even more importantly, this 
diagonalization significantly reduces the computational ef- 
fort required to solve the LU approximation and therefore 
transforms it into a more efficient method of numerically 
solving the three-dimensional Euler equations. 

1. Introduction 

Increasingly, attention is being focused upon implicit 
schemes as the need to develop more efficient numerical 
schemes becomes apparent. An implicit approximation is 
an attractive means of increasing the efficiency of the time- 
marching technique since an implicit scheme allows one to 
take much larger time steps than are normally allowed in 
most explicit methods. The drawback of most implicit 
schemes is that they are computationally more expensive, 
per time step, than explicit methods and thus interest is 
focused upon decreasing the computational effort required 
to solve the implicit approximation. Examples of schemes 
that address this issue are Chaussee and Pulliam’s Diago- 
nalized ADI 1 , Obayashi and Fujii’s LU-ADI 2 (based on the 
unsteady Euler and Navier-Stokes equations) and Yoon 
and Jameson’s LU-SSOR 3 (for the steady equations). 

The three-dimensional unsteady Euler equations can 
be discretized into an implicit approximation that, when 
written in delta-form, produces a large block banded ma- 
trix system that is impractical to solve without first be- 
ing approximately factored. In the present work, the LU 
Implicit Multigrid algorithm developed by Yokota and 
Caughey 4 for the three-dimensional unsteady Euler equa- 
tions has been made more efficient through a local diago- 
nalizing procedure that reduces the computational effort 
required to solve the LU approximation. This new Diag- 
onally Inverted LU scheme is related to the LU factoriza- 
tion first proposed by Jameson and Turkel 5 and applied 
to the two-dimensional Euler equations by Buratynski and 
Caughey 6 . The LU factorization produces two block trian- 
gular operators (one upper and one lower) which, through 
back substitution, can be solved by effectively explicit 
sweeps that require matrix inversions at every mesh cell in 
the domain. The efficiency of the LU scheme can be sig- 
nificantly increased by reducing the computational costs 


associated wdth these matrix inversions. This reduction 
can be achieved by a local diagonalizing transformation 
that decouples the matrix systems into scalar equations. 
Time conservation and stability are not altered by the de- 
coupling since the LU scheme’s differencing operators are 
not affected by the diagonalizing transformations. This re- 
sult is in contrast to the Diagonalized ADI method, which 
is produced at the expense of time accuracy. Thus the so- 
lution of the Diagonally Inverted LU scheme requires less 
computational effort than the original LU scheme and re- 
mains both unconditionally stable and time accurate. The 
Euler equations are numerically transformed to a gener- 
alized coordinate system and then discretized into a fi- 
nite volume approximation that is supplemented with an 
adaptive blend of second and fourth difference artificial 
dissipation. Numerical results show that the Diagonally 
Inverted LU scheme requires less than 1/4 the CPU time 
the original LU scheme would otherwise need. In addition, 
the diagonally Inverted LU scheme is developed within the 
framework of the multigrid method to accelerate the con- 
vergence rate for steady state calculations. Results for 
turbomachinery flows are presented to illustrate the Di- 
agonally Inverted LU scheme’s ability to calculate Euler-' 
flows and the convergence acceleration produced by the 
multigrid method. 

2. Analysis 


The unsteady three-dimensional Euler equations are 
written in divergence form and then transformed from the 
Cartesian coordinate system (x, y, z) to the generalized 
system (f,r?,£). The resulting equations can be written: 
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where p and P are the fluid density and pressure, (u, v, w ) 
are the Cartesian velocity components and E is tne total 
energy per unit volume. The total energy and pressure of 
a calorically perfect gas are related through the equation 
of state: 


p = (7 -,)( £ -^) 

The Jacobians of the coordinate transformation are as fol- 
lows: 
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where D is the determinant of the matrix J and the con- 
travariant velocities are defined as: 


face-centered central differences; J is the identity matrix; 
and A, B, and C are the flux Jacobian matrices relative 
to the vectors F, G, and H. 

3. Artificial Dissipation 

The finite volume formulation reduces to a central 
difference approximation on a uniform grid and therefore 
requires the addition of explicit artificial dissipation terms 
to suppress possible odd and even point oscillations and 
shock overshoots. Following the works of Jameson 7 and 
Pulliam 8 , fourth difference terms are added throughout 
the flow field to prevent odd-even decoupling while second 
difference terms are used to stablize the flow calculation 
near shocks. 

The conservatively added dissipative term is an adap- 
tive nonlinear blend of second and fourth differences that 
acts to turn on the second and turn off the fourth differ- 
ence terms near a shock. The dissipative terms added are 
of the form: 

fi = ~h« W Dh ~Wi + S t K (4) DS ((( jjWi (3) 

where a shock sensing term: 
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is used to define the nonlinear scaling factors: 
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The transformed Euler equations are discretized by 
a finite volume formulation that approximates the spatial 
differences as a net flux across the faces of each mesh cell. 
Global conservation and the admission of possible uniform 
flow solutions are insured by evaluating the flux vectors on 
the faces of the boundary-conforming mesh cells. This pro- 
cedure requires that both the flow variables and the flux- 
embedded geometric quantities be defined on the faces of 
the mesh cells during the flux evaluations although it is 
the cell-averaged flow variables that are calculated during 
the time and spatial marching. The geometric quantities 
are evaluated directly on the cell faces while the flow vari- 
ables are averaged over values found in adjacent cells. The 
unsteady equations can be discretized into an implicit ap- 
proximation that when written in a linearized delta form, 
produces a numerical scheme whose steady state solutions 
are independent of the time step size used in the time 
marching. The delta form, produced by linearizing the 
changes in the flux vectors through a Taylor series expan- 
sion about a time level n, can be written as follows: 

[/ + t(S ( A + + 6 ; C)\AW? ik = 

~At(^F + S n G + E,H)l. k ( 2 ) 

urn a|-a 

A W n = w n+ 1 -W n 


4 

K. tl = mai| 
*+ r 


( 0 \ 

V’ A >; + r i+} / 


where x ^ and x* 4 ' are scalar constants and At' is a one- 
dimensional unit Courant number time step scaling factor. 


4. LU Factorization 


The block-banded implict delta operator is factored 
into two block triangular operators to circumvent the large 
amounts of CPU time and temporary storage needed to 
solve the unfactored implicit approximation. The LU fac- 
torization, which is based on one-sided, implicit, spatial 
differences, can be written as follows: 
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where S + and 6~ are cell-centered forward and backward 
first differences and the flux Jacobian matrices are split 
and reconstructed as: 


and At is the time step size; 0 < p, < 1 is a parameter 
governing the degree of implicitness; S and 5 are cell- and 
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where 


5. Diagonal Inversion 


Ml = rriax(\X A \) 

Lb the maximum absolute-valued eigenvalue of the Jaco- 
bian matrix A, « 1 is a scalar constant governing the 
amount of implicit dissipation produced by the matix re- 
constructions and I is the identity matrix. This splitting 
is chosen to insure a nonsingular diagonal inversion and is 
similar to the Jameson and Turkel 6 splitting that produces 
diagonally dominant implicit factors for each matrix inver- 
sion in every mesh cell 4 . The implicit system of equations 
is solved in the following two steps: 

1) Lower Sweep 

[/ + Ai + 6^ Bi + 6~ Ci)j AY/j k = 

-A t(S t P + S n d + S f S + f) n iik ( 5 ) 

The right hand side of this system, the Residua], is evalu- , 
ated using the finite volume formulation of the numerical 
approximation and then inverted by an effectively explicit 
sweep through the domain in the positive (— ,r\— , di- 

rections. The resulting intermediate variables AY fn are 
then used to start the upper sweep. 

2) Upper Sweep 

[/ + pAt{6fA a + S+B 2 + 6?C 2 )] Atf,!}* = A Y£ k ( 6 ) 

An effectively explicit sweep through the domain is also 
required to solve this system, but in directions opposite 
to those taken in the lower sweep. The resulting flow field 
corrections are then used to update the flow field. 

W^ 1 = AW?- k + W? jk 

The LU factorization requires the solution of two 
block triangular operators each of which, through back 
substitution, can be reduced to simple 5x5 matrix sys- 
tems at every mesh cell. These reduced systems can be 
written: 

l) Lower Sweep 
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2) Upper Sweep 
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The solution of the LU factorization can be made 
more efficient by diagonally inverting its 5x5 matrix sys- 
tems. The diagonal inversion is based on the diagonalizing 
procedure used to classify the Euler equations formally as 
a hyperbolic system of partial differential equations. 

The Euler equations are classified as an hyperbolic 
system of partial differential equations since a similarity 
transformation exists such that: 

Q~ 1 (k 1 A + k 2 B + k 3 C)Q = k (9) 

where A is a diagonal matrix whose elements are real, and 
ki, k 2 , and k 3 are scalar constants. For the specific case of 
ki = k 2 = k 3 = 1, a corresponding diagonalizing transfor- 
mation can be derived from the general similarity trans- 
formation developed for the Euler equations in Cartesian 
coordinates by Warming, Beam, and Hyett 9 . The matrix 
Q has elements 
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and are typically inverted by Gaussian elimination. 
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and c is the local speed of sonnd. 

This diagonalizing transformation can be used to pro- 
duce the diagonal matrix 

Q~ 1 {A + B + C)Q = A 


which has elements 

A n = A 2 2 = A 3 3 = U + V + W 

x ii = u + v + w -cyjit + ii + q 

A n = U + V + W + cyjq + % + i% 

Aij = 0 when i ^ j 

Applying this local similarity transformation to the 
lower and upper sweeps produces the following scalar sys- 
tems: 
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2) Upper Sweep 
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The solution of these scalar equations requires sig- 
nificantly less computional effort than the original matrix 
equations. The diagonally inverted lower sweep requires 
approximately 855 operation counts per mesh cell as op- 
posed to the over 1453 operation counts per mesh cell 
required by the original lower sweep (a similar reduction 
is produced in the upper sweep). On a scalar machine, 
this reduction in operation counts translates into a 30-35 
percent reduction in CPU time, but on a vector machine 
the code can be restructured to produce an even greater 
increase in efficiency due to increased vectorization. All 
outward subroutine calls and short inner loops can be elim- 
inated from the lower and upper sweeps, thereby allow- 
ing most of the non-recursive parts of the code to vec- 
torize. This decrease in operation counts, together with 
the partial vectorization of the previously unvectorizable 
lower and upper sweeps, significantly increases the com- 
putational efficiency of the LU factorization. 


0. Initial and Boundary Conditions 
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The lower and upper sweeps have now been transformed 
into uncoupled systems of scalar equations that retain the 
stability and time accuracy of the original LU factoriza- 
tion. Unlike in the Diagonalized ADI method, time con- 
servation is not altered by the decoupling process since 
the similarity transformation is not factored out of the 
implicit spatial differences. For a time step defined as 


Initial conditions are needed to start the calculations 
and are chosen to be a uniform flow field based on an 
initial guess of the upstream flow. 

During the residual calculation, a no-flux condition 
is enforced at the solid boundaries by setting to zero the 
contravariant velocity component normal to the bound- 
ary. This treatment requires that only pressure needs to 
be specified along solid boundaries and is obtained from 
a three-dimensional interpretation of the normal momen- 
tum anaylsis developed by Rizzi 10 . For a solid boundary 
aligned with an r\ — f plane, the normal momentum anal- 
ysis requires that the following equation be satisfied. 
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where Cn is the Courant number used in the time march- 
ing, the vector components (m = 1, ...,5) of the lower and 
upper sweeps can now be written and solved as follows. 

1) Lower Sweep 

(Q- 1 Af^) m = 


where 

°i = (Dt x ) 2 + (L>£ y ) 2 + (L>6) 2 

0.2 = D 2 $ x r) x + D 2 ( y T]y + 

03 = D 2 Z X $ X + D^rjyty + I? 2 £*^ 

bl=pDV ^jm + jA£iA +w ?im^ 
h=pDW (jAm + „mA + jjk i) 

\ d* / 


4 



For a steady subsonic inflow boundary, total pressure, 
total temperature and two flow angles are specified while 
a one-dimensional Rieman invariant 


defined above requires 41/32 work units of computational 
effort where work units are normalized by a single flow 
calculation on the finest grid. 


R = q — 



8. Results 


is extrapolated from the interior flow field (similar to the 
treatment used by Chima 11 ). 

For steady subsonic outflow, the specified boundary 
condition is a nonreflective treatment, based on the work 
of Rudy and Strikwerda 12 , that attempts to minimizes un- 
wanted reflected waves from the outflow boundary. Static 
pressures are obtained from a radial momentum analysis 

1 d_P_ _ Ue_ 

p dr r 

(r is radial distance and Ujj is tangential velocity) and 
then coupled with the incoming compatibility relation to 
produce a nonreflective boundary condition 


dP 

dt 


pC dw 
dt 


+ a(P - Pref) = 0 
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where 


C = Cyjg + ^ 

P rc f is the static pressure obtained from the radial mo- 
mentum equation and a. e* 1 is a scalar constant. 

7. Steady State Calculations 

The calculation of steady state solutions is made more 
efficient by using local time-stepping and the multigrid 
method. 

Local time-stepping is used to optimize the time step 
throughout the flow field. A locally varing time step size, 
based on a constant Courant number, is used to create a 
warped time integration that can accelerate the calcula- 
tion to a steady state without affecting the steady state 
solution. 

The multigrid method is incorporated into the Diago- 
nally Inverted LU scheme to accelerate the removal of low 
frequency errors from the flow solution and thus increase 
the efficiency of the time-marching procedure. Following 
the work of Jameson 13 and Jameson and Yoon 14 , the flow 
solver is used to smooth out high frequency errors resolv- 
able on any current grid level, while the multigrid method 
is used to eliminate low frequency errors through a se- 
quence of flow calculations on coarser grids. The multigrid 
sequencing used is a three-level W- cycle that incorporates 
a double pass on the coarse grids. Coarse grid bound- 
ary conditions are identical to those used on the fine grid 
with the exception of the inflow/outflow conditions which 
are updated only on the fine grids. Coarse grid residu- 
als are kept smooth by adding only a constant coefficient 
second difference artificial dissipation term (the nonlin- 
ear blended terms are used only on the fine grid). This 
treatment attempts to limit the amount of high frequency 
errors reintroduced into the flow field by the upward inter- 
polation of the coarse grid corrections. The flow solver, in 
this case the Diagonally Inverted LU scheme, is invoked 
only once on each grid level and only before transfering 
the flow field to the next coarser grid. The multigrid cycle 


Numerical results are presented to illustrate the Di- 
agonally Inverted LU scheme’s ability to calculate three- 
dimensional turbomachinery flows and the convergence ac- 
celeration produced by the multigrid method. 

All results were calculated on H-type grids consist- 
ing of 64x16x16 mesh cells in the throughflow, blade-to- 
blade, and radial directions, respectively. The grids were 
generated using a modified version of the GRAPE code 15 
originally developed by Sorenson 16 . All calculations were 
performed on a CRAY X-MP, where a calculation consist- 
ing of 250 work units required approximately 0.9 million 
words of memory and 10 minutes of CPU time. This re- 
sult is in contrast to the original LU scheme for which a 
comparable calculation required over 42 minutes of CPU 
time. The diagonal inversion has significantly reduced the 
computational effort required to solve the LU factoriza- 
tion and has therefore transformed it into a more efficient 
numerical scheme. 

The first test case was a stationary (rot.= 0) cascade 
passage whose blade shapes were derived from the NACA 
0012 airfoil. The three-dimensional grid shown in Figure 
1 is based on a 36 blade axial rotor with a 0.825 hub- 
to-tip radius ratio and a 14.3 tip radius-to-blade chord 
ratio. Flow conditions of 0.66 exit static-to-inflow total 
pressure ratio (on the hub, P.R.= 0.66) and axial relative 
inflow (angy = 90, angz = 0) were used to produce a 
transonic flow through the turbomachinery passage. This 
case is used to illustrate the numerical scheme’s ability to 
calculate transonic flows and the convergence acceleration 
produced by the multigrid method. 

Convergence histories for calculations on a single grid 
and with the multigrid method are shown in Figures 2 
and 3. Both calculations were run at a Courant number 
of 10 and with local time-stepping. The logarithm of the 
flow field error and a measure of the developing super- 
sonic region are plotted versus the work units required to 
produce these results. The flow field error is character- 
ized by the logarithm of the rms rate of change of the 
continuity equation on the finest grid (scaled by its initial 
value), while the developing supersonic region (a measure 
of global convergence) is represented by the number of grid 
cells (NSUP) found to contain supersonic flow (scaled by 
its final value). The calculation performed on a single grid 
(Figure 2) shows a reduction in the average residual of al- 
most 3.5 orders of magnitude after 800 work units. The 
supersonic pocket has not yet converged to its final value 
after 800 work unit6, which suggests that the flow field 
error should be further reduced (Resl is the calculation’s 
starting residual while Res2 is the final residual). The 
multigrid method is incorporated into the time-stepping 
technique to produce a satisfactory reduction in the flow 
field error without having to resort to a large number of 
time steps. Figure 3 shows the convergence history of the 
three-level W-cycle, where a fully developed supersonic re- 
gion was produced in approximately 150 work units and 
4.5 order of magnitude drop in the flow field error was 
produced within 250 work units. The use of the multi- 
grid method has clearly accelerated the convergence rate 
of this steady state calculation. 

The blade surface and blade-to-blade Mach number 
distributions at mid-span are shown in Figures 4 and 5 
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to illustrate the accuracy and the shock-capturing capa- 
bilities of the Diagonally Inverted LU scheme. Near the 
blades the flow is accelerated to a peak Mach number of 
1.28 and then shocked down to subsonic flow through a 
shock that is captured over three mesh cells. 

The second test case is the Annular Cascade designed 
and extensively tested at NASA Lewis 17,18 . The computa- 
tional grid shown in Figure 6 is based on the full annular 
ring of 36 core turbine stator vanes. The geometry is a 
38.10 mm high untwisted blade of constant profile with 
an axial chord of 38.23 mm. The stator has a tip diameter 
of 508 mm and a 0.85 hub-to-tip radius ratio. 

Test conditions of ambient axial inflow and a 0.65 
hub-static to inlet-total pressure ratio produce a flow field 
with mean radius inlet and exit critical velocity ratios of 
0.231 and 0.778 respectively. To match the upstream flow 
conditions (an inflow Mach number of 0.211), the Euler 
calculations were run with a 0.66 hub-static to inlet-total 
pressure ratio. 

Figure 7 shows the convergence history of the Annular 
Cascade calculation (with multigrid) where a drop of 5 
orders of magnitude in the flow field error was produced 
within 250 work units. 

The resulting flow field is fully subsonic and compared 
with experimental data at three spanwise positions. Fig- 
ures 8, 9, and 10 compare the calculated blade surface 
static pressure distributions at 13.3, 50, and 86.7 percent 
span with the experimental data produced by Goldman 
and Seasholtz 18 . The trailing edge flow is not well cap- 
tured due to the bluntness of the blade’s trailing edge but 
the overall flow field agrees well with the experimental 
data. Figure 11 shows the blade-to-blade Mach number 
distribution at the mid-span position where the flow is 
turned approximately 67 degrees from the axial inflow. 
Figure 12 shows the critical velocity ratio distribution 
along the 50 percent radial span plane at the 150 percent 
axial chord location. The Euler equations cannot capture 
the viscous wake effects but the overall flow field does seem 
to be adequately modelled. The computational results 
compare well with the experimental data and demonstrate 
the scheme’s ability to calculate complex flows accurately. 

9. Concluding Remarks 

A new Diagonally Inverted LU Implicit scheme has 
been developed for the three-dimensional unsteady Euler 
equations. The matrix systems that are normally inverted 
in the LU scheme are treated by local diagonalizing trans- 
formations that decouple them, without a loss of time 
conservation, into systems of scalar equations. This de- 
coupling significantly reduces the computational effort re- 
quired to solve the LU approximation and therefore trans- 
forms it into a more efficient method of numerically solv- 
ing the unsteady three-dimensional Euler equations. Re- 
sults illustrate the scheme’s ability to calculate accurately 
three-dimensional flows through turbomachines and the 
convergence acceleration produced when coupled with the 
multigrid method. 
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Figure - 1. Computational grid for the NACA 0012 Cas- 
cade. 
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Figure - 3. Convergence history for the multigrid calcula- 
tion. 
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Figure - 2. Convergence history for the single grid calcu- 
lation. 


Figure - 4. Mach number blade distribution at 50 percent 
span. 
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Figure - 5. Mach number distribution at 50 percent span. 
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Figure - 7. Convergence history for the multigrid calcula- 
tion. 



Figure - 6. Computational grid for the Annular Cascade. 
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Figure - 8. Static pressure blade distribution at 13.3 per- 
cent span. 
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Figure - 11. Mach number distribution at 50 percent span. 


Figure - 9. Static pressure blade distribution at 50 percent 
span. 
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Figure - 10. Static pressure blade distribution at 86.7 per- 
cent span. 
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Summary 


A study of intelligent guidance and control concepts for protecting against the 
adverse effects of wind shear during aircraft takeoffs and landings has been begun. 
The principal objectives are to develop methods for assessing the likelihood of wind 
shear encounter, for deciding what flight path to pursue, and for using the aircraft’s 
full potential for combating wind shear. This study requires the definition of both 
deterministic and statistical techniques for fusing internal and external information, 
for making "go/no-go" decisions, and for generating commands to the aircraft’s 
cockpit displays and autopilot for both manually controlled and automatic flight. 

The program has begun with the development of a real-time expert system for 
pilot aiding that is based on the results of the FAA Windshear Training Aids Pro- 
gram. A two-volume manual that presents an overview, pilot guide, training pro- 
gram, and substantiating data provides guidelines for this initial development. The 
Expert System to Avoid Wind Shear (ESAWS) currently contains over 140 rales and 
is coded in the LISP programming language for implementation on a Symbolics 
3670 LISP Machine. 

This six -month progress report includes a brief introduction to ESAWS, as 
well as numerous appendices that describe the logic and code of the program in con- 
siderable detail. This development is especially fluid, and these materials are pre- 
sented as evidence of progress made to date rather than as documentation of a 
working system; consequently, the logic is subject to further change, and no guaran- 
tee is made that the submitted coding is "bug-free" or even representative of 
ESAWS's final structure. Program logic is being developed by Professor Stengel 
and graduate student Alexander Stratton; Mr. Stratton has been responsible for all 
program coding and has authored a considerable portion of this progress report. 
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1. LISP Implementation of Wind Shear Expert System 

Wind Shear Expert System software has been written in Common LISP for a 
Symbolics 3670 LISP machine, using elementary list operations. Execution of the 
system and input/output operations are carried out in the standard LISP environ- 
ment, using a "LISP Listener". Currently, the user acts as both flight crew and 
simulation environment; additions planned for the system, including a FORTRAN 
aircraft simulation, will add realism for system development and demonstration and 
will complete the model of flight crew decision-making and control presented in the 
FAA Windshear Training Aid. 

Software for the Wind Shear Expert System consists of parameter definitions 
(Appendix A), rule bases (Appendix B), and functional LISP procedures 
(Appendices C to F). Rule bases are defined as LISP variables using the function 
"defvar " . A rule base is simply a list of related rules (see Appendix B for English 
translations of rules). In the code, each mle is a list containing seven elements. The 
first element of a rule is the mle's name. Two other elements, called the"if-part" or 
premise and the "then-part" or consequent , are executable LISP statements that are 
evaluated if the mle is tried. Two of the elements are character strings that give an 
English translation of the mle; these strings are typed to a LISP Listener if the 
premise is evaluated and does not return "NIL". 

The two remaining elements of a rule are lists of parameters. One list consists 
of the variables whose values are set if the consequent is evaluated. This list is a 
"flag" for a goal-directed search procedure. The other parameter list consists of 
variables that must have values for the premise to be evaluated; this list is a "flag" 
for a forward-chaining procedure and is not presently used. Editing features of the 
Symbolics machine enable nested lists such as the mle bases to be easily entered into 
files in a neat format that simplifies identification and later modification. Grouping 
the rules into separate mle bases breaks down the expert system functions into sub- 
functions in a natural manner. 

Expert system parameters store data for the expert system, and are objects of 
search procedures (see Appendix A). Parameters have been defined using 
" defvar ", and are bound using " setq " in mles and functional procedures. Some 
parameters are additionally given property values, to expand them to vectors. 

Some parameters have been grouped to aid variable binding. With hindsight, this 
has worked well for internal parameters that are given values during search, but it 
is cumbersome for representing other parameters external to the system. 
Redefinition of external data, perhaps using "defflavor" and "make-instance", may 
have advantages and is being considered. 


Functional LISP procedures retrieve rule elements, evaluate rules, conduct 
goal-directed searches, and bind variables. They are defined using "defun" and are 
called from a LISP Listener or from within the rule bases. The goal-directed search 
procedure " get-value- of " is currently the only search procedure being used, 
although other procedures have been written. "Get-value-of\ is called with a goal 
parameter and a rule base as arguments. It searches the rule base for rules that have 
a "flag" indicating they can set the value of the goal parameter and tries them 
usingthe procedure "tryrule". "Get-value-of is called to begin execution of the 
expert system and is called frequently within the rule base, resulting in backward- 
chaining. The development of LISP procedures is essentially complete; however, 
FORTRAN procedures for simulating, estimating, and controlling aircraft and sys- 
tems must be implemented. The Symbolics LISP machine permits interaction of 
LISP and FORTRAN procedures, making the demonstration of FA A Windshear 
Training Aid procedures and other algorithms possible. 

The expert system may be called by a LISP Listener or by a FORTRAN 
program; the former method has been used up until now. The expert system 
presently sends documentation to a LISP Listener while it is operating. It prompts 
the user for inputs through a LISP Listener; the user can give additional inputs by 
temporarily halting execution of the expert system. This interim arrangement has 
been adequate for development thus far, but it does not adequately demonstrate the 
interaction of the expert system with its operating environment. 

To adequately demonstrate the system and aid further development, a simula- 
tion environment must be implemented that includes aircraft, on-board systems, 
ground-based information sources, and "mother nature". Modelling events outside 
of the system adds autonomy and realism to the simulation. To better understand the 
interaction of the expert system with its environment, a multi-window display must 
be developed that gives information about several aspects of the system and its envi- 
ronment simultaneously. The Symbolics machine has programs that are intended to 
streamline the development of multi- window displays; our intention is to quickly 
develop a simple means of presenting a variety of information. 


2. Search Procedures for Wind Shear Expert System 

The Wind Shear Expert System, in its current state, is principally a model of 
pilot and crew decision-making and control, as described in the FAA Windshear 
Training Aid. The four primary functions of the system are: 

1. MONITORING — Observe sensors, receive reports, alerts, warnings 

2. ASSESSMENT ~ Detect wind shear encounters, determine if there 

are signs of wind shear, and if it is safe to continue 

3. PLANNING — Determine what actions and precautions are taken 

4. ACTION — In automatic mode, execute standard, recovery, and 

go-around procedures; in semi-automatic mode, issue 
commands to flight directors 

These functions are performed as side effects of a goal-directed search for para- 
meter values in a set of rules. A mle, in the context of this report, is an IF-THEN 
statement; a list of the 141 mles implemented on a Symbolics 3670 LISP machine 
is enclosed (see Appendix B). A search procedure is simply a process of "trying" 
rules in a selective fashion. "Trying" a rule means that the premise of a rule is 
examined to determine its truth or falsehood, and the consequent is examined if the 
premise is true. Often the consequent is a direction to perform some action or pro- 
cedure. Typically, the value of a parameter is set in the consequent. The search 
procedure is given a goal parameter, it finds mles that have consequents that set the 
value of the goal parameter, and it tries them sequentially, stopping when a mle is 
tried whose premise is true. If another parameter, whose value has not been deter- 
mined, is encountered in a rule's premise, that parameter becomes the new goal 
parameter and a new search commences for the value of that parameter. If none of 
the mles capable of setting a goal parameter have premises that are true, the user is 
asked to give a value for the parameter in question. This enables the expert system 
to be tested with an incomplete mle base. 

To demonstrate the actions of the search procedure, consider the following 
example. The expert system is invoked periodically from a simulation by directing 
the search procedure to get a value for the parameter " RBC -search-complete " . 
Turning to the section "Parameters for Wind Shear Expert System", we fmd that 
"RBC -search-complete" is set in the EXECUTIVE mle base. Now turning to 
Appendix B, we find two mles that can set this parameter: 


"IF the mission phase is not landed, 

THEN we must complete executive procedures, 
and the rule-based controller search is not complete. 

IF the mission phase is landed, 

THEN the rule-based controller search is complete." 

Trying the first rule, we must consider the truth or falsehood of the expression, "the 
mission phase is not landed". To know this, we must find out what the mission phase 
is. Turning to the list of parameters, we find a parameter " mission-phase " that is set 
in the MISSION PHASE rule base. The rules within the MISSION PHASE rule base 
are tried next. When the mission phase has been determined, we return to consider- 
ing the truth or falsehood of the first premise. If, for example, we determined that 
the mission phase was not " landed ", we would then consider the consequent of the 
first rule, "we must complete executive procedures, and the rule-based controller 
search is not complete." The first part of the consequent is actually a call to the 
search procedure (rules were translated into English to facilitate readability; regret- 
tably, this can make it more difficult to follow the search process), to get a value for 
the parameter "executive-search-complete". Once this search is complete, the 
parameter "RBC-search-complete" is set to NIL and the top-level search is finished. 


3. Plans for Improvement of Wind Shear Expert System 

The Wind Shear Expert System has developed primarily from the top level 
downward, branching out as the number of rules increases. This approach, one 
of the advantages of the expert system paradigm for control system development, 
tends to leave open-ended parameters on the lowest level. The material enclosed 
does not represent a complete rule-based controller and should be viewed as an 
interim development. Currently, there are a few parameters that do not have 
any rules or procedures to give them values. Additionally, there are procedures 
called by the expert system that are not yet implemented as software. Briefly, a list 
of major items still needing attention includes: 

• Simulation environment 

• Improved user interface 

• Improved relevance information and testing 

• Control system architecture 

• Target trajectory generation 

• Prediction procedures for aircraft and weather 

• Rules to interface with new procedures 

Addition of the above items should complete an expert system model of the FAA 
Windshear Training Aid and provide the basis for further development of a rule- 
based pilot aid for wind shear survival. 


Appendix A 

Parameters for Wind Shear Expert System 

Below is a list of the current parameters contained in the wind shear expert 
system. Parameters have been classified according to how they are given values. 
These classifications are: 

INTERNAL (I) - Parameter is internal to the expert system and gets its 

value from a goal-directed search procedure. 

PRESET (Pr) - Parameter is set by flight crew, through the pilot 

interface. 

STATE VARIABLE (S) - Parameter is a state estimate, set by a state estimator. 

OUTSIDE (O) - Parameter's value is given to the program by an 

outside interface. 

1. Internal Parameters 

PARAMETER NAME TYPE VALUES (INTERPRETED) 


EXECUTIVE RULEBASE sets: 

RBC-search-complete I T (we should keep searching) 

NIL (we should stop searching) 

executive-procedures -complete I T (procedures are not complete) 


MISSION PHASE RULEBASE sets: 


mission-phase 


current-airport 

are-we-committed 

WINDSHEAR RULEBASE sets: 
microburst-encounter 

DETECTION RULEBASE sets: 
detection-complete 


I PREFLIGHT-DECISION stage 
TAKEOFF (pilot committed) 
CLIMBOUT (off the ground) 
CRUISE (takeoff completed) 
PRELANDING-DECISION stage 
APPROACH (pilot committed) 
FLARE (below flare height) 
LANDING (on runway) 
LANDED (operations completed) 

Pr/I [any airport name] (nearest 
airport) 

Pr/I T (flight crew committed to 
takeoff/Iand) 

NIL (flight crew not committed) 


I T (we are having a microburst 
encounter) 

NIL (no microburst encounter) 


I T (search of detection rule base 
completed) 


DEVIATION RULEBASE sets: 


aircraft-response I [list] contains one or more of: 

MICROBURST (significant 
deviations): 

ABOVE-TARGET-AIRSPEED 
BELOW-TARGET-AIRSPEED 
ABOVE-TARGET- VERTICAL- 
SPEED 

BELOW-TARGET- VERTICAL- 
SPEED 

ABOVE-TARGET-PITCH- 

ATTITUDE 

BELOW-TARGET-PITCH- 

ATTITUDE 

ABOVE-TARGET- 

GLIDESLOPE-DEV. 

BELOW-TARGET- 

GLIDESLOPE-DEV. 

ABOVE-TARGET-THROTTLE- 

POS. 

BELOW-TARGET-THROTTLE- 

POS. 


target-airspeed I 

airspeed-deviation I 

target-vertical-speed I 

vertical-speed-deviation I 

target-pitch-attitude I 

attitude-deviation I 

target-glideslope-displacement I 


[numerical] (knots) 

[numerical] (deviation from target, 
knots) 

[numerical] (ft/min) 

[numerical] (deviation from target, 
ft/min) 

[numerical] (deg) 

[numerical] (deviation from target, 
deg) 

[numerical] (dots) 


• 

glideslope-dev iation 

I 

[numerical] (deviation from target, 
dots) 


throttle-deviation 

I 

[numerical] (deviation from target, 

%) 


COMMUNICATION RULEBASE sets: 


new-information-received 

I 

T (new information is received) 
NIL (no new info received) 


incident-reported 

I 

T (if there was an incident, it’s 
reported) 


tower-informed-goa 

I 

T (tower informed of intent to go- 
around) 

NIL (tower not yet informed of 
intent) 


tower-informed-delay 

I 

T (tower informed of intent to 
delay) 

NIL (tower not yet informed) 


ACTION RULEBASE sets: 




actions-taken 


I T (actions have been taken) 


STANDARD PROCEDURES 

RULEBASE sets: 


standard-procedures-complete 

I 

T (standard procedures are now 
complete) 


configuration-set-for-takeoff 

I 

T (configuration is set for takeoff) 
NIL (configuration not set for 
takeoff) 


configuration-set-for-landing 

I 

T (configuration is set for landing) 
NIL (configuration is not set for 
landing) 
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glideslope-established Pr/I T (glideslope has been established) 

NIL (glideslope is not yet 
established) 

RECOVERY PROCEDURES RULEBASE sets: 

recovery-procedures -complete I T (recovery procedures are now 

complete) 

GO-AROUND PROCEDURES RULEBASE sets: 

go-around-procedures-complete I T (go-around procedures are 

complete) 

DELAY PROCEDURES RULEBASE sets: 

delay-procedures -complete I T (delay procedures are complete) 

PLANNING RULEBASE sets: 

recommended-procedures I T (mission plan is now updated) 

RUNWAY SELECTION RULEBASE sets: 

recommended-runway I [runway id] (new runway 

• selection) 

SPEED SELECTION RULEBASE sets: 

• recommended- Vr I [numerical] (precautionary 

rotation speed) 

recommended-Va I [numerical] (precautionary 

approach speed) 
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FLAP SELECTION RULEBASE sets: 


recommended-flaps I [numerical] (precautionary flap 

setting) 

ASSESSMENT RULEBASE sets: 

recommended-action I CONTINUE operations as planned 

CON S3DER-PREC AUTION S 
DELAY-ALTER-ABORT 
EXECUTE-RECOVERY 
procedures 


RISK RULEBASE sets: 

risk I [numerical] (overall risk factor for 

wind shear - values correspond to: 

0 - no probability of wind shear 

1 - low probability of wind shear 

2 - medium probability of shear 

3 or greater - high probability 


WEATHER RISK RULEBASE sets: 

weather I [numerical] (risk of wind shear 

derived from weather reports, e.g., 
SIGMETs, same severity scale as 
"risk") 

suspected-runways I [list of runway id's] (runways 

reported or observed to have 
conditions suggesting wind shear) 


LOCALIZED WINDS RULEBASE sets: 


localized-strong- winds I [numerical] values correspond to: 

0 - no indications of wind shear 

1 - indications, off the flight path 

2 - indications, convectirig onto 

path 

3 - indications of wind shear on 

path 


HEAVY PRECIPITATION RULEBASE sets: 

heavy-precipitation I [numerical] (same scale as 

localized-strong winds) 

RAIN SHOWER RULEBASE sets: 

rainshower I [numerical] (same scale as 

localized-strong winds) 


LIGHTNING RULEBASE sets: 

lightning I [numerical] (same scale as 

localized-strong winds) 

VIRGA RULEBASE sets: 

virga I [numerical] (same scale as 

localized-strong winds) 


TURBULENCE RULEBASE sets: 

turbulence I [numerical] (moderate or greater 

same scale as localized-strong- 
winds) 


2. Preset Parameters 


origin-airport 

Pr 

[any airport name] (origin airport) 

taxi-speed-limit 

Pr [numerical] (taxi speed limit, 
knots) 

destination-airport 

Pr 

[any airport name] (next 
destination) 

flare-agl 

Pr 

[numerical] (flare altitude, ft) 

Vr 

Pr 

[numerical] (standard rotation 
speed for this configuration, 
knots) 

Va 

Pr 

[numerical] (standard approach 
speed for this configuration, 
knots) 

critical-agl 

Pr 

[numerical] (height below which 
ground contact is major concern, 
ft) 

aircraft-type 

Pr 

[any type in PWG] (type of 
aircraft) 

VI 

Pr 

[numerical] (abort groundspeed, 
knots) 

suitable-runways 

Pr 

[list of runways] (suitable 
runways) 

V r-f ield-length-limit 

Pr 

[numerical] (field length limit 
rotation speed for this 
configuration, knots) 


S [numerical] (alt. above ground, ft) 


3. State Variables 

agl 


airspeed 

vertical-speed 

pitch-attitude 

glideslope-displacement 

groundspeed 

runway-remaining 

4. Outside Parameters 

wind shear-alert 


traffic 

pirep* 


llwas* 


forecast* 

temp/dewpt-spread 


S 

S 

s 

s 

s 

s 


o 

o 

o 

o 

o 

o 


[numerical] (current airspeed, ft) 

[numerical] (ft/min) 

[numerical] (deg) 

[numerical] (dots) 

[numerical] (groundspeed, knots) 

[numerical] (current est. of 
available runway, ft) 


T (on-board wind shear detection 
sys. alert) 

NIL (no on-board wind shear det. 
sys. alert) 

T (threatening traffic ahead) 

NIL (no threatening traffic ahead) 

[numerical] (amount of speed 
lost/gained by some other flight 
crew, knots) 

[numerical] (amount of speed 
differential that triggered some 
llwas alert, knots) 

[list] (forecast of convective 
weather) 

[numerical] (degrees F) 
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ATIS* 


ASWW* 


SIGMET* 


onboard-radar* 


tower-report* 


TDWR* 


wind-profiler* 


0 TURBULENCE noted in report 
RAINSHOWERS noted in report 
LIGHTNING noted in report 
VIRGA noted in report 
HEAVY-PRECIPITATION noted 
LLWS (low-level wind shear 
noted) 

O TURBULENCE noted in report 
RAINSHOWERS noted in report 
LIGHTNING noted in report 
VIRGA noted in report 
HEAVY-PRECIPITATION noted 
LLWS (low-level wind shear 
noted) 

O TURBULENCE noted in report 
RAINSHOWERS noted in report 
LIGHTNING noted in report 
VIRGA noted in report 
HEAVY-PRECIPITATION noted 
LLWS (low-level wind shear 
noted) 

O TURBULENCE observed 
HEAVY-PRECIPITATION 
observed 

O TURBULENCE noted by tower 
RAINSHOWERS noted by tower 
LIGHTNING noted by tower 
VIRGA noted by tower 
HEAVY-PRECIPITATION noted 
WINDSHEAR noted by tower 

O TURBULENCE noted by TDWR 
HEAVY-PRECIPITATION noted 
WINDSHEAR noted by TDWR 

O TURBULENCE detected by 
profiler 

WINDSHEAR detected by profiler 


forward-looking * 

0 

TURBULENCE detected 
RAINSHOWERS detected 
LIGHTNING detected 
VIRGA detected 
HEAVY-PRECIPITATION 
detected 

WINDSHEAR detected 

observation 

0 

[numerical] (phenomenon seen by 
flight crew on flight path, values 
correspond to: 

0 - no suspicious phenomenon seen 

1 - localized strong winds seen 

2 - heavy precipitation seen 

3 - rain showers seen 

4 - lightning seen 

5 - virga seen 

LSS* 

0 

LIGHTNING detected by lightning 
sensor 


* The parameter also has properties indicating airports and runways to which it 
pertains. 


Appendix B 

Rule Bases for Wind Shear Expert System 


EXECUTIVE RULEBASE 

IF the mission phase is not landed, 

THEN we must complete executive procedures, 

and the rule-based controller search is not complete. 

IF the mission phase is landed, 

THEN the rule-based controller search is complete. 

IF the mission phase is preflight decision stage or prelanding decision stage, 
THEN we must determine recommended procedures, 
and executive procedures are complete. 

IF the mission phase is takeoff, 

or the mission phase is climbout, 
or the mission phase is approach, 
or the mission phase is flare, 
or the mission phase is landing rollout, 

THEN we must take control actions, 

and executive procedures are complete. 

IF the mission phase is takeoff complete, 

THEN the current airport is now the destination airport, 
and the pilot is not committed to land, 
and executive procedures are complete. 


MISSION PHASE RULEBASE 

IF the current airport is the origin airport, 
and our altitude(ft) is 0 feet, 
and the pilot is not committed to takeoff, 
THEN the mission phase is preflight decision stage. 


IF the current airport is the origin airport, 
and our altitude (ft) is 0 feet, 
and the pilot is committed to takeoff, 

THEN the mission phase is takeoff. 

IF the current airport is the origin airport, 
and our altitude (ft) is 0 feet, 

and our airspeed (knots) is greater than the taxi speed limit (knots), 
and the flight crew has not committed to takeoff, 

THEN the flight crew has passively committed to takeoff, 
and the mission phase is takeoff. 

IF the current airport is the origin airport, 

and our altitude (ft) is greater than 0 feet, 
and our altitude (ft) is less than 1000 feet, 

THEN the mission phase is climbout. 

IF the current airport is the origin airport, 

and our altitude (ft) is greater than 1000 feet, 

THEN the mission phase is takeoff complete. 

IF the current airport is the destination airport, 

and our altitude (ft) is greater than 1000 feet, 
and the flight crew is not committed to land, 

THEN the mission phase is prelanding decision stage. 

IF the current airport is the destination airport, 
and our altitude (ft) is below 1000 feet., 
and we are not flight crew is not committed to land, 

THEN the flight crew has passively committed to land, 
and the mission phase is approach. 

IF the current airport is the destination airport, 
and our altitude (ft) is less than 1000 feet, 
and our altitude (ft) is above flare altitude (ft), 
and the flight crew is committed to land, 

THEN the mission phase is approach. 

IF the current airport is the destination airport, 

and our altitude (ft) is below flare altitude (ft), 
and the flight crew is committed to land, 

THEN the mission phase is flare. 
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IF the current airport is the destination airport, 
and our altitude (ft) is 0 feet, 

and our airspeed (knots) is above the taxi speed limit (knots), 
THEN the mission phase is landing roll. 

IF the current airport is the destination airport, 
and our altitude (ft) is 0 feet, 

and our airspeed (knots) is below the taxi speed limit (knots), 
THEN the mission phase is landed. 


WINDSHEAR RULEBASE 

IF the search of the detection rulebase is complete, 

and the aircraft trajectory is deviating significantly from the target 
trajectory, 

THEN there is a microburst encounter. 

IF a wind shear alert is sounded by an on-board wind shear detection system, 

THEN there is a microburst encounter. 

IF a microburst encounter has already begun, 

and the aircraft trajectory is deviating somewhat from the target trajectory, 

THEN there is a microburst encounter. 

IF a microburst encounter has not begun, 

and the aircraft trajectory is not deviating significantly from the target 
trajectory, 

THEN there is no microburst encounter. 


DETECTION RULEBASE 

IF the mission phase is approach or prelanding decision stage, 
THEN airspeed deviations must be monitored, 

and vertical speed deviations must be monitored, 
and pitch attitude deviations must be monitored, 
and glideslope displacement deviations must be monitored, 
and throttle position deviations must be monitored, 
and the search of the detection rulebase is complete. 


IF the mission phase is neither approach nor prelanding decision stage, 

THEN airspeed deviations must be monitored, 

and vertical speed deviations must be monitored, 
and pitch attitude deviations must be monitored 
and the search of the detection rulebase is complete. 

DEVIATION RULEBASE 

IF airspeed (knots) is below target airspeed (knots) by more than 15 knots, 

THEN the aircraft is deviating significantly from the target trajectory. 

IF airspeed (knots) is above target airspeed (knots) by more than 15 knots, 

THEN the aircraft is deviating significantly from the target trajectory. 

IF airspeed (knots) is below target airspeed (knots) by between 5 and 15 knots, 

THEN the aircraft is deviating somewhat from the target trajectory, 
and the type of deviation includes "below target airspeed". 

IF airspeed (knots) is above target airspeed (knots) by between 5 and 15 knots, 

THEN the aircraft is deviating somewhat from the target trajectory, 
and the type of deviation includes "above target airspeed". 

IF vertical speed (ft/min) is below target vertical speed (ft/min) by more than 500 
ft/min, 

THEN the aircraft is deviating significantly from the target trajectory. 

IF vertical speed (ft/min) is above target vertical speed (ft/min) by more than 500 
ft. /min., 

THEN the aircraft is deviating significantly from the target trajectory. 

IF vertical speed (ft/min) is below target vertical speed (ft/min) by between 200 and 
500 ft./min., 

THEN the aircraft is deviating somewhat from the target trajectory, 
and the type of deviation includes "below target vertical speed". 

IF vertical speed (ft/min) is above target vertical speed (ft/min) by between 200 and 
500 ft./min., 

THEN the aircraft is deviating somewhat from the target trajectory, 
and the type of deviation includes "above target vertical speed". 

IF pitch attitude (deg) is below target pitch attitude (deg) by more than 5 degrees, 

THEN the aircraft is deviating significantly from the target trajectory. 


IF pitch attitude (deg) is above target pitch attitude (deg) by more than 5 degrees, 

THEN the aircraft is deviating significantly from the target trajectory. 

IF pitch attitude (deg) is below target pitch attitude (deg) by between 2 and 5 
degrees, 

THEN the aircraft is deviating somewhat from the target trajectory, 
and the type of deviation includes "below target pitch attitude". 

IF pitch attitude (deg) is above target pitch attitude (deg) by between 2 and 5 
degrees, 

THEN the aircraft is deviating somewhat from the target trajectory, 
and the type of deviation includes "above target pitch attitude". 

IF glideslope displacement (dots) is below target glideslope displacement (dots) by 
more than 1 dot, 

THEN the aircraft is deviating significantly from the target trajectory. 

IF glideslope displacement (dots) is above target glideslope displacement (dots) by 
more than 1 dot, 

THEN the aircraft is deviating significantly from the target trajectory. 

IF glideslope displacement (dots) is below target glideslope displacement (dots) by 
between 0.4 and 1 dot, 

THEN the aircraft is deviating somewhat from the target trajectory, 

and the type of deviation includes "below target glideslope displacement". 

IF glideslope displacement (dots) is above target glideslope displacement (dots) by 
between 0.4 and 1 dot, 

THEN the aircraft is deviating somewhat from the target trajectory, 

and the type of deviation includes "below target glideslope displacement". 


COMMUNICATION RULEBASE 

IF the outside interface indicates new information, 

THEN new information has been received. 

IF the outside interface does not indicate new information, 
THEN no new information has been received. 


IF there was a microburst encounter, 

and there is no longer a microburst encounter, 
THEN the incident must be recorded now. 

IF the above premise is not true, 

THEN there is no need to report an incident. 


ACTION RULEBASE 

IF we are having a microburst encounter, 
and guidance mode is automatic, 

THEN recovery procedures must be completed, 
and control actions are taken. 

IF we are not having a microburst encounter, 
and guidance mode is automatic, 
and recommended procedures have been determined, 
and standard procedures are recommended, 

THEN standard procedures must be completed, 
and control actions are taken. 

IF we are not having a microburst encounter, 
and guidance mode is automatic, 
and go-around procedures are recommended, 

THEN go-around procedures must be completed, 
and control actions are taken. 

IF we are not having a microburst encounter, 

and delaying procedures are recommended, 

THEN delaying procedures must be completed, 
and control actions are taken. 

IF we are having a microburst encounter, 
and guidance mode is semi-automatic, 

THEN semi-automatic recovery procedures must be completed, 
and control actions are taken. 


IF we are not having a microburst encounter, 
and guidance mode is semi-automatic, 
and recommended procedures have been determined, 
and standard procedures are recommended, 

THEN semi-automatic standard procedures must be completed, 
and control actions are taken. 

IF we are not having a microburst encounter, 
and guidance mode is semi-automatic, 
and go-around procedures are recommended, 

THEN semi-automatic go-around procedures must be completed, 
and control actions are taken. 


STANDARD PROCEDURES RULEBASE 

IF the mission phase is takeoff, 

and the configuration has not been set for the takeoff roll, 
THEN the configuration must be set for the takeoff roll, 
and standard procedures are complete. 

IF the mission phase is takeoff, 

and the airspeed (knots) is below Vr (knots), 

THEN takeoff control laws must be executed, 
and standard procedures are complete. 

IF the mission phase is takeoff, 

and the airspeed (knots) is above Vr (knots), 

THEN rotation control laws must be executed, 
and standard procedures are complete. 

IF the mission phase is climbout, 

THEN climbout control laws must be executed, 
and standard procedures are complete. 

IF the mission phase is approach, 

and the configuration has not been set for landing, 

THEN the configuration must be set for landing, 
and standard procedures are complete. 


IF the mission phase is approach, 

and the glideslope has not been established, 

THEN we must complete procedures to establish the glideslope, 
and standard procedures are complete. 

IF the mission phase is approach, 

and the glideslope has been established, 

THEN approach control laws must be executed, 
and standard procedures are complete. 

IF the mission phase is flare, 

THEN flare control laws must be executed, 
and standard procedures are complete. 

IF the mission phase is landing roll, 

THEN landing control laws must be executed, 
and standard procedures are complete. 


RECOVERY PROCEDURES RULEBASE 

IF the altitude (ft) is above critical (ft), 
and the aircraft type is an L-1011, 
and the mission phase is takeoff, 

THEN controls must be configured for a recovery, 

and the desired thrust setting is maximum rated, 
and the desired pitch attitude (deg) is 17.5 degrees, 
and recovery procedures are complete. 

IF the altitude is above critical, 

and the aircraft is not an L-1011, 

or the aircraft is an L-1011 and the mission phase is approach, 
THEN controls must be configured for a recovery, 
and the desired thrust setting is maximum rated, 
and the desired pitch attitude (deg) is 15 degrees, 
and recovery procedures are complete. 


IF the altitude is below critical, 

and the mission phase is not takeoff or landed, 

THEN controls must be configured for a recovery, 
and the desired thrust setting is overboost, 
and the desired pitch attitude (deg) is 20 degrees, 
and recovery procedures are complete. 

IF the mission phase is takeoff, 

and the groundspeed (knots) is less than VI (knots), 
THEN controls must be configured for takeoff abort, 
and the desired thrust setting is reverse, 
and recovery procedures are complete. 

IF the mission phase is takeoff, 

and the groundspeed (knots) is greater than VI (knots), 
and the airspeed (knots) is less than Vr (knots), 
and the runway remaining (ft) is greater than 2000 feet, 
THEN controls must be configured for takeoff, 

and the desired thrust setting is maximum rated, 
and recovery procedures are complete. 

IF the mission phase is takeoff, 

and the airspeed (knots) is greater than Vr (knots), 
THEN controls must be configured for a recovery, 

and the desired thrust setting is maximum rated, 
and the desired pitch attitude (deg) is 15 degrees, 
and recovery procedures are complete. 

IF the mission phase is takeoff, 

and the runway remaining (ft) is less than 2000 feet, 
THEN controls must be configured for a recovery, 
and the desired thrust setting is maximum rated, 
and the desired pitch attitude (deg) is 15 degrees, 
and recovery procedures are complete. 


GO-AROUND PROCEDURES RULEBASE 


IF the tower has not been informed of our intent to go-around, 
and there is no traffic in our way, 

THEN go-around control laws must be executed, 

and the tower must be informed of our intent to go-around, 
and go-around procedures are complete. 

IF the tower has not been informed of our intent to go-around, 
and there is traffic in our way, 
and our altitude is above critical, 

THEN approach control laws must be executed, 

and the tower must be informed of our intent to go-around, 
and go-around procedures are complete. 

IF the tower has not been informed of our intent to go-around, 
and there is traffic in our way, 
and our altitude is below critical, 

THEN go-around control laws must be executed, 

and the tower must be informed of our intent to go-around, 
and go-around procedures are complete. 

IF the tower has been informed of our intent to go-around, 

THEN go-around control laws are in effect, 

and go-around procedures are complete. 


DELAY PROCEDURES RULEBASE 

IF the tower has not been informed of our intent to delay, 

and our ETA at the airport (min) is less than 20 minutes, 
THEN the tower must be informed of our intent to delay, 
and delay procedures are complete. 

IF the tower has been informed of our intent to delay, 

THEN delay procedures are complete. 


PLANNING RULEBASE 


IF the recommended action is to delay, alter, or abort, 
and the mission phase is preflight decision stage, 
or the mission phase is prelanding decision stage, 

THEN delaying procedures are recommended. 

IF the recommended action is to delay, alter, or abort, 
and the mission phase is takeoff, 
or the mission phase is climbout, 
or the mission phase is flare, 
or the mission phase is landing roll, 

THEN standard procedures are recommended. 

IF the recommended action is to delay, alter, or abort, 
and the mission phase is approach, 
and the altitude (ft) is above decision height (ft), 

THEN go-around procedures are recommended. 

IF the recommended action is to delay, alter, or abort, 
and the mission phase is approach, 
and the altitude (ft) is below decision height (ft), 

THEN standard procedures are recommended. 

IF the recommended action is to take precautions, 

and precautions have already been taken, or the flight crew has rejected 
precautions, 

THEN standard procedures are recommended. 

IF the recommended action is to take precautions, 

and the mission phase is preflight decision stage, 

THEN we must determine the best runway, 

and we must determine the best flap setting (deg), 

and we must determine the best Vr (knots), 

and the precautions must be accepted or rejected by the flight crew. 

IF the recommended action is to take precautions, 

and the mission phase is prelanding decision stage, 

THEN we must stabilize approach at 1000 feet, 

and we must determine the best flap setting (deg), 

and we must determine the best Va (knots), 

and the precautions must be accepted or rejected by the flight crew. 
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IF the recommended action is to take precautions, 

and the mission phase is not preflight or prelanding decision stage, 

THEN standard procedures are recommended. 

IF the recommended action is to continue, 
and the mission phase is approach, 
and altitude (ft) is below critical altitude (ft), 
and the glideslope has not been established, 

THEN go-around procedures are recommended. 

IF the recommended action is to continue, 
and the mission phase is approach, 
and the altitude (ft) is above the decision altitude (ft), 
and the weighted norm of the 2-sigma touchdown dispersion (ft) is greater 
than the dispersion tolerance (ft), 

THEN go-around procedures are recommended. 

IF the recommended action is to continue, 
and the above premises are not true, 

THEN standard procedures are recommended. 


RUNWAY RULEBASE 

IF some of the suitable runways are suspected of wind shear, 

THEN the best runway is the longest of the suitable runways not suspected. 

IF none of the suitable runways are suspected of wind shear, 

THEN the best runway is the longest of the suitable runways. 


SPEED RULEBASE 

IF field-length-limit Vr (knots) is greater than the standard Vr (knots), 
and is less than 20 knots greater than the standard Vr (knots), 

THEN the best Vr (knots) is the field-length-limit Vr (knots). 

IF field-length-limit Vr (knots) is more than 20 knots greater than the standard Vr 
(knots), 

THEN the best Vr (knots) is 20 knots plus the standard Vr (knots). 


IF this rule is fired, 

THEN the best Va (knots) is 20 knots plus the standard Va (knots). 


FLAP RULEBASE 

IF the mission phase is preflight decision stage, 
and the aircraft type is a B727, 

THEN the best flap setting (deg) is 15 degrees. 

(72 more rules like the above, for different AlC configurations ) 

ASSESSMENT RULEBASE 

IF there is no new information received, 

or new information has been received, 

and all of the mles in the risk rulebase have been tried, 

and the overall risk factor is 0 (none), 

THEN the recommended action is to continue. 

IF the overall risk factor is 3 or greater (high), 

THEN the recommended action is to delay, alter, or abort. 

IF the overall risk factor is 2 (medium), 

THEN the recommended action is to consider precautions. 

IF the overall risk factor is 1 (low), 

THEN the recommended action is to continue. 


RISK RULEBASE 

IF this rule is fired, 

THEN the overall risk factor is set to 0 
and the list of causes is set to nil. 

IF there is an indication of convective weather on our flight path, 
THEN all of the rules in the weather risk rulebase must be tried, ' 
and the weather risk factor is added to the overall risk factor, 
and weather is added to the list of causes. 


IF there is a pilot report of speed loss (knots) less than 15 knots (but nonzero), 
and this report is relevant to our current flight path, 

THEN the overall risk factor is increased by 2, 

and the current runway is added to the list of suspected runways, 
and the pirep is added to the list of causes. 

IF there is a pilot report of speed loss (knots) greater than 15 knots, 
and this report is relevant to our current flight path, 

THEN the overall risk factor is increased by 3, 

and the current runway is added to the list of suspected runways, 
and the pirep is added to the list of causes. 

IF there is a LLWAS alert with a speed differential (knots) less than 20 knots, 
and this alert is relevant to our current flight path, 

THEN the overall risk factor is increased by 2, 

and the current runway is added to the list of suspected runways, 
and the LLWAS alert is added to the list of causes. 

IF there is a LLWAS alert with a speed differential (knots) greater than 20 knots, 
and this alert is relevant to our current flight path, 

THEN the overall risk factor is increased by 3, 

and the current runway is added to the list of suspected runways, 
and the LLWAS alert is added to the list of causes. 

IF there is a pilot report of speed loss (knots) less than 15 knots (but nonzero), 
and this report is near but not on our current flight path, 

THEN the overall risk factor is increased by 1 , 

and the current runway is added to the list of suspected runways, 
and the pirep is added to the list of causes. 

IF there is a pilot report of speed loss (knots) greater than 15 knots, 
and this report is near but not on our current flight path, 

THEN the overall risk factor is increased by 2, 

and the current runway is added to the list of suspected runways, 
and the pirep is added to the list of causes. 

IF there is a LLWAS alert with a speed differential (knots) less than 20 knots, 
and this alert is near but not on our current flight path, 

THEN the overall risk factor is increased by 1, 

and the current runway is added to the list of suspected runways, 
and the LLWAS alert is added to the list of causes. 


IF there is a LLWAS alert with a speed differential (knots) greater than 20 knots, 
and this alert is near but not on our current flight path, 

THEN the overall risk factor is increased by 2, 

and the current runway is added to the list of suspected runways, 
and the LLWAS alert is added to the list of causes. 

IF there is a forecast of convective weather on our flight path, 

THEN the overall risk factor is increased by 1, 
and forecast is added to the list of causes. 


WEATHER RISK RULEBASE 
IF this rule is fired, 

THEN the weather risk factor is reset to zero. 

IF there are indications of localized strong winds on our flight path, 

THEN the weather risk factor is increased by 3, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of heavy precipitation on our flight path, 

THEN the weather risk factor is increased by 3, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of a rain shower on our flight path, 

THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of lightning on our flight path, 

THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of virga on our flight path, 

THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of moderate or greater turbulence on our flight path, 
THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 


IF the ATIS report indicates a temperature/dew point spread (° F) greater than 
30 OF, 

THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of localized strong winds near but not on our flight path, 

THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of heavy precipitation near but not on our flight path, 

THEN the weather risk factor is increased by 2, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of rain showers near but not on our flight path, 

THEN the weather risk factor is increased by 1, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of lightning near but not on our flight path, 

THEN the weather risk factor is increased by 1, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of virga near but not on our flight path, 

THEN the weather risk factor is increased by 1, 

and the affected runways are added to the list of suspected runways. 

IF there are indications of moderate or greater turbulence near our flight path, 

THEN the weather risk factor is increased by 1, 

and the affected runways are added to the list of suspected runways. 


LOCALIZED WINDS RULEBASE 

IF ATIS information indicates low-level wind shear, 
and the wind shear is on our flight path, 
or an Aviation Severe Weather Watch indicates low-level wind shear, 
and the wind shear is on our flight path, 
or a SIGMET indicates low-level wind shear, 
and the wind shear is on our flight path, 
or a tower report indicates low-level wind shear, 
and the wind shear is on our flight path, 
or a Terminal Doppler Weather Radar indicates low-level wind shear, 
and the wind shear is on our flight path, 


or a wind profiler indicates low-level wind shear, 
and the wind shear is on our flight path, 
or a forward-looking system indicates low-level wind shear, 
and the wind shear is on our flight path, 
or a flight crew observation indicates low-level wind shear on our flight path, 
THEN there are indications of localized strong winds on our flight path. 

IF ATIS information indicates low-level wind shear, 

and the wind shear is near but not on on our flight path, 
or an Aviation Severe Weather Watch indicates low-level wind shear, 
and the wind shear is near but not on our flight path, 
or a SIGMET indicates low-level wind shear, 

and the wind shear is near but not on our flight path, 
or a tower report indicates low-level wind shear, 

and the wind shear is near but not on our flight path, 
or a Terminal Doppler Weather Radar indicates low-level wind shear, 
and the wind shear is near but not on our flight path, 
or a wind profiler indicates low-level wind shear, 

and the wind shear is near but not on our flight path, 
or a forward-looking system indicates low-level wind shear, 
and the wind shear is near but not on our flight path, 
or a flight crew observation indicates low-level wind shear near our flight path, 
THEN there are indications of localized strong winds near our flight path. 


HEAVY PRECIPITATION RULEBASE 

IF ATIS information indicates heavy precipitation, 

and the heavy precipitation is on our flight path, 
or an Aviation Severe Weather Watch indicates heavy precipitation, 
and the heavy precipitation is on our flight path, 
or a SIGMET indicates heavy precipitation, 

and the heavy precipitation is on our flight path, 
or the on-board radar indicates heavy precipitation, 
and the heavy precipitation is on our flight path, 
or a tower report indicates heavy precipitation, 
and the heavy precipitation is on our flight path, 
or a Terminal Doppler Weather Radar indicates heavy precipitation, 
and the heavy precipitation is on our flight path, 
or a forward-looking system indicates heavy precipitation, 
and the heavy precipitation is on our flight path, 
or a flight crew observation indicates heavy precipitation on our flight path, 
THEN there are indications of heavy precipitation on our flight path. 


IF ATIS information indicates heavy precipitation, 

and the heavy precipitation is near but not on on our flight path, 
or an Aviation Severe Weather Watch indicates heavy precipitation, 
and the heavy precipitation is near but not on our flight path, 
or a SIGMET indicates heavy precipitation, 

and the heavy precipitation is near but not on our flight path, 
or the on-board radar indicates heavy precipitation, 

and the heavy precipitation is near but not on our flight path, 
or a tower report indicates heavy precipitation, 

and the heavy precipitation is near but not on our flight path, 
or a Terminal Doppler Weather Radar indicates heavy precipitation, 
and the heavy precipitation is near but not on our flight path, 
or a forward-looking system indicates heavy precipitation, 

and the heavy precipitation is near but not on our flight path, 
or a flight crew observation indicates heavy precipitation near our flight path, 
THEN there are indications of heavy precipitation near our flight path. 


RAINSHOWER RULEBASE 

IF ATIS information indicates a rain shower, 
and the rain shower is on our flight path, 
or an Aviation Severe Weather Watch indicates a rain shower, 
and the rain shower is on our flight path, 
or a SIGMET indicates a rain shower, 

and the rain shower is on our flight path, 
or the on-board radar indicates a rain shower, 
and the rain shower is on our flight path, 
or a tower report indicates a rain shower, 
and the rain shower is on our flight path, 
or a Terminal Doppler Weather Radar indicates a rain shower, 
and the rain shower is on our flight path, 
or a forward-looking system indicates a rain shower, 
and the rain shower is on our flight path, 
or a flight crew observation indicates a rain shower on our flight path, 

THEN there are indications of a rain shower on our flight path. 

IF ATIS information indicates a rain shower, 

and the rain shower is near but not on on our flight path, 
or an Aviation Severe Weather Watch indicates a rain shower, 
and the rain shower is near but not on our flight path, 
or a SIGMET indicates a rain shower, 


and the rain shower is near but not on our flight path, 
or the on-board radar indicates a rain shower, 

and the rain shower is near but not on our flight path, 
or a tower report indicates a rain shower, 

and the rain shower is near but not on our flight path, 
or a Terminal Doppler Weather Radar indicates a rain shower, 
and the rain shower is near but not on our flight path, 
or a forward-looking system indicates a rain shower, 
and the rain shower is near but not on our flight path, 
or a flight crew observation indicates a rain shower near our flight path, 
THEN there are indications of a rain shower near our flight path. 


LIGHTNING RULEBASE 

IF ATIS information indicates lightning, 

and the lightning is on our flight path, 
or an Aviation Severe Weather Watch indicates lightning, 
and the lightning is on our flight path, 
or a SIGMET indicates lightning, 

and the lightning is on our flight path, 
or the Lightning Sensor System indicates lightning, 
and the lightning is on our flight path, 
or a tower report indicates lightning, 
and the lightning is on our flight path, 
or a forward-looking system indicates lightning, 
and the lightning is on our flight path, 
or a flight crew observation indicates lightning on our flight path, 

THEN there are indications of lightning on our flight path. 

IF ATIS information indicates lightning, 

and the lightning is near but not on on our flight path, 
or an Aviation Severe Weather Watch indicates lightning, 
and the lightning is near but not on our flight path, 
or a SIGMET indicates lightning, 

and the lightning is near but not on our flight path, 
or the Lightning Sensor System indicates lightning, 
and the lightning is near but not on our flight path, 
or a tower report indicates lightning, 

and the lightning is near but not on our flight path, 
or a forward-looking system indicates lightning, 
and the lightning is near but not on our flight path, 
or a flight crew observation indicates lightning near our flight path. 


THEN there are indications of lightning near our flight path. 


VIRGA RULEBASE 

IF ATIS information indicates virga, 

and the virga is on our flight path, 
or an Aviation Severe Weather Watch indicates virga, 
and the virga is on our flight path, 
or a SIGMET indicates virga, 

and the virga is on our flight path, 
or a tower report indicates virga, 
and the virga is on our flight path, 
or a forward-looking system indicates virga, 
and the virga is on our flight path, 
or a flight crew observation indicates virga on our flight path, 
THEN there are indications of virga on our flight path. 

IF ATIS information indicates virga, 

and the virga is near but not on on our flight path, 
or an Aviation Severe Weather Watch indicates virga, 
and the virga is near but not on our flight path, 
or a SIGMET indicates virga, 

and the virga is near but not on our flight path, 
or a tower report indicates virga, 

and the virga is near but not on our flight path, 
or a forward-looking system indicates virga, 

and the virga is near but not on our flight path, 
or a flight crew observation indicates virga near our flight path, 
THEN there are indications of virga near our flight path. 


TURBULENCE RULEBASE 

IF ATIS information indicates moderate or greater turbulence, 
and the turbulence is on our flight path, 

or an Aviation Severe Weather Watch indicates moderate or greater turbulence, 
and the turbulence is on our flight path, 
or a SIGMET indicates moderate or greater turbulence, 
and the turbulence is on our flight path, 
or the on-board radar indicates turbulence, 
and the turbulence is on our flight path, 
or a tower report indicates turbulence. 


and the turbulence is on our flight path, 

or a Terminal Doppler Weather Radar indicates moderate or greater turbulence, 
and the turbulence is on our flight path, 
or a wind profiler indicates moderate or greater turbulence, 
and the turbulence is on our flight path, 
or a forward-looking system indicates moderate or greater turbulence, 
and the turbulence is on our flight path, 
or a flight crew observation indicates moderate or greater turbulence on our 
flight path, 

THEN there are indications of moderate or greater turbulence on our flight path. 

IF ATIS information indicates moderate or greater turbulence, 
and the turbulence is near but not on on our flight path, 
or an Aviation Severe Weather Watch indicates moderate or greater turbulence, 
and the turbulence is near but not on our flight path, 
or a SIGMET indicates moderate or greater turbulence, 
and the turbulence is near but not on our flight path, 
or the on-board radar indicates turbulence, 

and the turbulence is near but not on our flight path, 
or a tower report indicates turbulence, 

and the turbulence is near but not on our flight path, 
or a Terminal Doppler Weather Radar indicates moderate or greater turbulence, 
and the turbulence is near but not on our flight path, 
or a wind profiler indicates moderate or greater turbulence, 
and the turbulence is near but not on our flight path, 
or a forward-looking system indicates moderate or greater turbulence, 
and the turbulence is near but not on our flight path, 
or a flight crew observation indicates turbulence near our flight path, 

THEN there are indications of moderate or greater turbulence near our flight path. 


Appendix C 

LISP Code for Rule Base Executive Functions 
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variable TOWER-REPORT has properties winds 


Variable ONBOARD-RADAR can have values NIL, TURBULENCE, and HEAVY-PRECIPITATION 
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11 and this is a strong indication of a microburst encounter, 

(setq aircraft-response '(microburst)) ' air speed -deviation) 
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(recvryl " ACTION: In microburst, AGL is above critical, aircraft type is L-1011, 


APPLY PULL THRUST IMMEDIATELY 
PULL TO 17.5 DEGREES PITCH ATTITUDE 
RESPECT STICK SHAKER 
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AIRSPEED IS ABOVE VR 
APPLY FULL THRUST IMMEDIATELY 
PULL TO 15 DEGREES PITCH ATTITUDE 
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( delay 1 " ACTION: Decision stage, tower has not been informed of our intention to delay, 

(eq tower-informed-delay nil) 'tower-informed-delay 


" so the tower will be informed of our intent to delay takeoff." 

(get-value-of 1 tower- informed-delay comm-rules) 'delay-procedures-complete) 
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LISP Code for Rules of Wind Shear Assessment 
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setq weather-rules 
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E-3 


( weather9 " ASSESSMENT: There is a relevant tower report of heavy precipitation 

(and (eq tower-report 'blowing-dust) 

(neq ( is-it-relevant 'tower-report current-runway) nil)) 'tower-report 
" and this is a significant indication of a microburst." 

(and (setq weather (+ 3 weather)) 
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INTELLIGENT 
GUIDANCE AND CONTROL 
FOR WINDSHEAR ENCOUNTERS 


INTRODUCTION 

AND 

PROGRESS REPORT 


Alex Stratton 
Princeton University 
June 23-24, 1988 


Overview 


• Introduction 


Expert System Development 


Plans for Future Work 


Motivation for the Study 


Strong windshears, especially microbursts, can cause 
fatal accidents 


Pilot proficiency in coping with microburst windshear is 
difficult given rarity of phenomenon 


Proper decision making and control strategy can enhance 
the possibility of avoidance and survival 


F.A.A. Pilot Windshear Guide establishes such a strategy 


Intelligent G&C Research Goals 


Implement the Pilot Windshear Guide pilot and 
crew model as a rule-based (intelligent) controller 


Investigate guidance and control concepts applied 
to windshear encounters 


Develop a rule-based pilot aid for windshear survival 


Rule-based Pilot Aid - Schematic 




EXTERNAL 


Weather reports 
LLWAS 
PIREPs 
SIGMETS 
Tower reports 

NEXRAD 

TDWR 


LANGUAGE 

INTERFACE 


MESSAGE . 
INTER- 
PRETATION 
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Sources of Decision and Control 

Knowledge 


Prior knowledge of windfield - Optimal control laws 


No knowledge of windfield - Strategy based on current 
knowledge in cockpit - Pilot Windshear Guide 


Pilot Windshear Guide 


Completed as part of the FAA's Windshear Training 
Aid 


Reference guide for flight crews detailing windshear 
avoidance, recognition, and recovery techniques 


Provides a model of pilot and crew decision-making 
that increases flight safety 




I 

Rule-based Systems for Control 

» 


Federated systems - Symbolic supervisory control 

Separate numeric processing 



Integrated systems - Combined symbolic and numeric control 


INTEGRATED SYSTEM j 

■>!vS‘X< , >>>> V ' VA ^ V '^' AW ^ A ^ Ww ^ % ^‘^ V ^ AVjW ‘ V VVWAVyW.VVA*A'A%\VW//AVAV.V^\VfcV/A*Xv'. , '-‘*! , /:v.:v'.'’' 



Combined 


Combined 

11 


Symbolic 

xSS :• iv.’V: : 

Symbolic 



and 

:ftS¥K»!5A?SSte 

and 

'■ 


Numerical 

, :•: • ' 
..^ : 

Numerical 

i-iS-l 

:•• .- : - 
\ •. : 

Processing 


Processing 

y : : : : 

f | 



IpIlPlIti 
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Expert System Control Functions 


• MONITORING 

• ASSESSMENT - INTERPRETATION - DIAGNOSIS 

• PREDICTION - PLANNING 


• ACTION - CONTROL 
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Monitoring and Assessment 


MONITORING- 

Observe onboard sensors 

Receive reports, alerts, warnings 

Send reports, communicate with flight crew 


ASSESSMENT - 
Determine the risk of continuing 
Decide if the flight plan should be altered 
Decide if windshear has been encountered 


> 


> 


Planning and Action 


• PLANNING - 

Delay or abort terminal operations if necessary 
Take additional safety precautions 


# ACTION - 

Select and implement guidance and control laws 
Send commands to control surfaces and throttle 
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Expert System Implementation 


• LISP — Knowledge representation, symbolic processing 


• FORTRAN — Numerical procedures, algorithms 


• Symbolics 3670 LISP Machine 


Flow of Control Driven By 
Goal-Directed Search Procedure 


Basic search procedure evaluates rules to set the value 
of a goal parameter 


Search procedure is often involked in the premise of 
a rule, resulting in backward-chaining 


More complicated search procedures are built with the 
basic procedure 


Graphical Representation of Search 
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Example of Knowledge Acquisition 

PILOT WINDSHEAR GUIDE STATEMENT: 

"The choice of takeofff flap selection is dependant 
on the airplane type. The following should be con- 
sidered unless limited by obstacle clearance and/or 
climb gradient: 

Airplane Takeoff 

Type Flap Setting 

B727 15 

B737 5 to 15 


PSEUDO-LISP TRANSLATION: 

IF recommended-action equals "take-precautions" 
AND aircraft-type equals "B727" 

THEN set recommended- flap-setting to "15" 


H-17 


Expert System Detailed Structure 










Next Additions to Expert System 


FORTRAN procedures and simulation environment 
will be completed next 


Low-level rules to interface FORTRAN procedures will 
complete knowledge base for Pilot Windshear Guide 


Improved user interface will aid the understanding and 
further development of expert system pilot aid 


Next Additions to Expert System 


\ 
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Plans for Future Work 


• Consider impact of next-generation sensors and ATC 


• Real-time implementation of RBC 


• Consider pilot/RBC interface 
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